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Preface 


Preface 

The  purpose  of  this  study  was  to  develop  a  technique  to  linearize  a 
non-linear  turbofan  engine  into  linear  state  space  format.  The  need  for 
this  study  arose  from  the  desire  to  better  understand  state-of-the-art 
engine  control  by  development  of  accurate  linear  state  space  models. 

To  develop  the  linear  state  space  models.  I  used  a  non-linear 
turbofan  engine  simulation,  commonly  referred  to  as  a  cycle  deck.  The 
turbofan  engine  simulation  used  in  this  study  was  for  the  F 1 0 1  engine  in 
use  on  the  B-1B  aircraft.  Although  the  F 1 0 1  engine  is  used  in  this 
particular  study,  the  methodology  for  the  development  of  linear  state 
space  models  is  valid  for  other  turbofan  engines. 

I  would  very  much  like  to  thank  Dr.  Robert  Calico  for  his  patience 

and  guidance,  not  just  on  the  writing  of  this  thesis,  but  for  all  of  the 

classes  I  have  taken  from  him  and  all  the  knowledge  and  insights  he  has 
given  me.  I  have  drawn  on  that  knowledge  many  times  and  will  continue 
in  the  future.  Thank  you.  I  would  also  like  to  thank  my  branch, 

ASD/ENFPA,  in  allowing  me  the  use  of  computer  resources  and  time  to 

complete  this  study.  In  particular,  Russel  Denney.  His  keen 
understanding  of  engine  simulation  models  was  a  constant  source  for  me 
in  understanding  the  intricate  coding  that  makes  up  a  non-linear  engine 


simulation  model . 


Gregory  L.  Thelen 
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Abstract 

Present  generation  turbofan  engines  use  hydro/mechanical  control 
governors  to  regulate  fuel  flow  and  control  engine  performance. 
State-of-the  -art  and  future  engines  will  have  the  ability  to  use 
digital  control  logic  to  control  engine  performance.  Due  to  these 
advances  in  engine  control  capability,  there  is  a  need  to  linearly  model 
the  turbofan  engines  and  develop  control  systems  to  understand  and 
optimize  engine  performance.  This  paper  will  describe  the  means  of 
developing  linear  state  space  models  which  model  transient  turbofan 
engine  performance. 

The  F101  turbofan  engine,  used  on  the  B-1B  bomber,  will  be  used  as 
the  example  with  the  linear  state  space  models  being  derived  from  the 
non-linear  F101  engine  computer  simulation  model.  The  internal 
convergence  logic  of  the  F101  engine  simulation  will  be  used  to  derive 
the  individual  elements  making  up  the  linear  state  space  models.  The 
linear  state  space  models  will  consist  of  both  high  speed  and  low  speed 
rotor  dynamics  and  turbine  inlet  temperature  heat  soak  dynamics.  State 
space  Inputs  considered  will  be  fuel  flow  and  engine  exit  nozzle  area. 
Also  dicussed  in  this  paper  will  be  linear  analyitlc  equations  in  3tate 
space  format  and  their  comparative  accuracies  to  the  models  derived 
using  the  F101  non-linear  computer  simulation  model. 

Based  on  the  linear  state  space  models  developed  in  this  paper, 
control  systems  will  be  designed  and  implemented  into  the  F101  engine 
computer  model.  Transient  performance  will  be  compared  between  current 
engine  control  design  and  the  control  design  based  on  the  linear  state 


LINEAR  STATE  SPACE  MODELING  OF  A  TURBOFAN  ENGINE 

I .  Introduction 

The  dynamic  response  of  turbofan  engine  systems  may  be  accurately 
analyzed  with  modern  state  space  techniques.  Recent  advances  in  engine 
technology,  such  as  the  extensive  use  of  variable  geometry  and  digital 
electronic  controls,  are  ideally  suited  to  state  space  analysis. 
However,  developement  of  a  state  space  model  can  be  challanging.  The 
difficulty  in  developing  linear  models  lie3  in  the  use  of  non-linear 
thermodynamic  computer  engine  models,  commonly  referred  to  as  cycle 
decks.  Cycle  decks  are  usually  purchased  by  the  Air  Force  from  the 
engine  manufacturer.  Engine  modeling  techniques  can  vary  greatly  from 
manufacturer  to  manufacturer  which  makes  it  very  difficult  to  develop  a 
generic  method  of  using  the  non-linear  cycle  decks  from  the  different 
engine  manufacturers  in  deriving  linear  state  space  models.  There  are, 
however,  some  common  denominators  in  engine  thermodynamic  modeling  among 
manufacturers  that  serve  as  a  basis  in  the  development  of  linear  state 
space  models  from  non-linear  cycle  decks. 

This  paper  will  first  examine  how  non-linear  cvcle  decks 
thermodynamically  model  engine  dynamics  using  a  simple  turbojet  as  an 
example.  The  results  will  be  analogous  to  the  more  complex  turbofan.  A 


non-linear  turbofan 

cvcle  deck 

will  then  be 

used 

to  derive  linear 

models . 

Also  discussed  in  this 

paper  will  be 

the 

use 

of  linear 

analytic 

equat 1 ons 

to  model  the  thermod vnami c 

proc 

messes 

within  the 

engine.  Finally,  using  the  developed  linear  models,  a  rate  feedback 
design  which  takes  advantage  of  present  digital  technology  will  be  shown 
to  improve  engine  response  characteristics  over  present  day 
hydro/mechanical  engine  control  designs  used  in  the  F101/B-1B  turbofan 


II  Non- 1 1  near  Cycle  Deck 


A  non-linear  cycle  deck  thermodynamically  balances  pressures, 
temperatures,  and  flows  across  various  engine  components,  the  most 
common  being  the  engine  inlet,  low  pressure  fan  (LPF)  ,  high  pressure 
compressor  (HPC) ,  high  pressure  turbine  (HPT),  the  low  pressure  turbine 
(LPT)  and  the  exhaust  nozzle.  The  cycle  deck  satisfies  both  energy  and 
flow  continuity  requirements  based  on  inlet,  fan,  compressor,  combustor, 
turbine  and  nozzle  component  character istics  predefined  from 
experimental  testing.  The  component  characteristics  are  mapped 
(defined)  as  functions  of  temperature,  pressure,  and  mass  flow  rates  and 
integrated  into  the  thermodynamic  processes  and  calculations  internal  to 
the  cycle  deck  computor  simulation.  Therefore,  the  cycle  deck  engine 
simulation  models  individual  components  with  unique  temperature, 
pressure  and  flow  characteristics  and  through  energy  and  flow  continuity 
relationships  establish  their  interdependence.  To  simplify  the 
discussion  of  internal  thermodynamic  balancing  of  a  non-linear  engine 
cycle  deck,  a  simple  turbojet  will  be  considered  . 

Turbojet  Engine 

Figure  1  is  a  representation  of  a  turbojet  engine  consisting  of  a 
compressor,  combustor,  turbine  and  nozzle.  All  pressures  and 

temperatures  are  absolute  unless  otherwise  indicated,  and  the  inlet 
considered  lsentropical ly  ideal.  Also  shown  are  component 

characteristic  mappings  of  the  compressor,  turbine.  and  nozzle. 
Combustor  efficiency  (P^/P  )  is  a  function  of  temperature,  pressure,  and 
geometry  of  the  combustor  and  is  considered  defined  from  a  table  look-up 


once  the  thermodynamic  states  are  known  at  Station  3. 

The  power. lever  angle  (PLA) ,  or  engine  throttle  position,  defines 

the  demanded,  corrected,  engine  shaft  speed,  Nfc.  In  the  F101/B-1B  and 

the  F110/F-15  turbofan  engines,  the  schedule  is  imposed  on  a  mechanical 

cam  linked  to  the  PLA  position.  For  a  given  PLA  setting,  the  engine 

will  operate  somewhere  along  a  constant  corrected  speed  line,  N^,  at 

some  operating  point  defined  by  the  compressor  map  (see  Figure  1).  Any 

point  along  a  given  will  define  a  mass  flow  rate,  pressure  ratio,  and 

efficiency  across  the  compressor.  An  initial  guess  is  made  by  the 

engine  cycle  deck  program  as  to  where  along  the  N^line  the  engine  will 

operate  defining  the  corresponding  pressure  ratio,  P^/P^,  corrected 

mass  flow  rate,  m„,_,  ,  defined  as  m„_  VT  -  .  /P„_ .  and  efficiency,  r)  ,  across 
/OK  is  is  is  c 

the  compressor.  It  is  assumed  that  inlet  conditions,  T„_  and  Pn_  are 

25  l  5 

always  known  from  instrumentation  in  the  inlet.  An  expression  for  the 
compressor  work  in  terms  of  temperature  differences  across  the 
compressor  and  based  on  the  guessed  operating  point  on  the  compressor 
map  is  given  in  Equation  (1)  as, 


r  -  i 


T  -  T 
25  3 


where  ^  is  a  function  of  temperature  and  fuel/air  ratio./.  Compressor 
work  and  are  now  defined. 

Thermodynamic  properties  equating  T^and  T^across  the  combustor  can 
be  written  as, 


(1  *  /,Cp4T4  =  Cp3T3  +  'V 
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VVV V  /  V 


where,  Cp  -  Specific  heat 

r)  =  Combustor  efficiency 

Q  =  Heating  value  of  fuel. 

The  cycle  deck  guesses  a  value  for  /  and  given  T3  from  Equation  (1), 

Equation  (2)  may  be  solved  for  T^.  Recalling  the  values  for  ,  defined 

by  PLA  position,  and  m^,  defined  from  /  and  m^ ,  P^,  defined  from  P^  and 

rv  ,  and  T,,  the  turbine  pressure  ratio,  P./P-,  can  be  determined  from 
b  4  4  o 

the  turbine  map  characteristics  (Figure  1).  Turbine  work  can  now  be 
calculated  using  the  equation  below. 


T4  '  T5  =  V4  1  - 


P  /  P 

4  5 


For  steady  state  operation,  the  work  being  done  across  the 
compressor  must  be  equal  to  the  work  across  the  turbine  or  acceleration 
in  rotor  speed  will  occur.  If,  based  on  the  guessed  fuel/air  ratio,  /, 
the  compressor  and  turbine  work  are  not  equal,  fuel/air  ratio  will  be 
used  to  iterate  between  Equations  (2)  and  (3)  until  the  the  compressor 
and  turbine  work  are  balanced. 

The  nozzle  oressure  ratio,  P_/P  ,  and  the  nozzle  corrected  mass  flow 

5  o 

rate.  m. VtT/P- .  can  be  calculated  from  the  following  equations: 


PK/P  *  (P./P  )(P,/PJ(P./P,)IP./P.) 
5o  2o  32  43  54 


where  Pq  is  the  ambient  pressure  ,  and 


"V  T  5 


m.  t  T 
4 


p 

4  r  4 
P 

*  c  ^ 
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All  the  ratios  on  the  right  hand  side  of  the  above  equations  are 

known.  However,  P./P  and  the  corrected  mass  air  flow  must  also  lie 

o  o 

somewhere  on  the  nozzle  map  (Figure  1).  Continuity  must  be  satisfied 

between  calculated  corrected  mass  flow  rate  and  nozzle  corrected  mass 

flow  rate  characteristics.  Recall  that  calculated  values  for  P_/P  and 

5  O 

corrected  mass  flow  rate  are  based  on  the  initially  guessed  compressor 
operating  point.  If  the  calculated  nozzle  corrected  mass  flow  rate  and 
the  corrected  mass  flow  rate  based  on  the  nozzle  map  characteristics  do 
not  match  for  a  given  P_/P  ,  a  new  compressor  operatng  point  along  the 

speed  line  must  be  chosen  until  continuity  is  satisfied.  Once 

satified.  the  entire  thermodynamic  cycle  is  satisfied  for  both  energy 
and  continuity. 

Various  techniques  for  iteration  on  fuel/air  ratio  and  compressor 
operating  line  exist  and  vary  from  simplified  techniques  which  can  be 
accomplished  on  a  hand-held  calculator  to  more  accurate  methods  which 
require  the  Newton-Raphson  technique  such  as  is  currently  used  in 
General  Electric  engine  cycle  simulations*,  but  the  basic  premise  of 
energy  and  continuity  satisfaction  must  always  be  met. 


Turbofan  Engine 


Consider  a  turbofan  as  a  high  rotor  speed  turbojet  surrounded  by 
another  low  rotor  speed  turbojet.  Stage  by  stage  the  thermodynamic 
energy  equations  are  satisfied  in  the  same  manner  as  the  simple  turbojet 
example.  Figure  2  shows  a  typical  turbofan  engine.  The  additional 
components  distinguishing  it  from  a  turbojet  are  the  low  pressure  fan, 
(LPF) ,  low  pressure  turbine,  (LPT),  and  the  bypass  duct.  Like  the 
turbojet,  inlet  conditions  are  known.  Work  across  the  LPF  and  LPT  must 


be  equal  just  as  work  across  the  high  pressure  compressor,  (HPC) , 


T  % 


•X 


V, 

f 


■ri 


and 


the  high  pressure  turbine,  (HPT).  Also  like  the  turbojet,  an  operating 
point  for  tie  LPF  is  guessed  by  assuming  a  bypass  ratio  (ratio  of  air 
through  the  duct  divided  by  the  air  through  the  compressor)  and  low 
rotor  speed.  This  serves  to  define  a  pressure  ratio  and  corrected  mass 
air  flow  across  the  LPF  and  in  turn,  inlet  conditions  into  the  HPC. 
Continuity  is  satisfied  by  requiring  static  pressures  of  the  duct  and 
high  speed  ror,or  mass  flow  rates  at  the  mixing  point  of  the  duct  and 
core  be  equal.  The  mass  air  flows  mix  shortly  after  station  5.  Static 
pressures  will  be  equal  provided  no  shocks  exist  at  the  mixing  point 
which  is  normally  the  case.  The  F101/B-1B  assures  this  to  be  the  case 
by  scheduling  duct  Mach  number  as  a  function  of  engine  inlet 
temperature.  Duct  mass  flow  velocity  is  controlled  by  varying  nozzle 
exit  area.  Nozzle  flow  continuity  of  the  combined  duct  and  core  mass 
air  flow  must  also  be  satisfied  based  on  nozzle  map  characteristics. 

The  energy  and  continuity  conditions  for  a  turbofan  are  satisfied  in 
an  analagous  manner  to  those  of  a  turbojet.  The  only  real  conceptual 
difference  is  in  the  number  of  energy  and  continuity  conditions  to  be 
satisfied.  The  important  point  to  realize  is  that  in  equilibrium,  flow 
continuity  is  satisfied  along  with  the  work  energy  conditions.  Any 
torque  applied  off  equilibrium  (off-torque)  on  either  rotor  would  force 
a  rematching  of  energy  (spool  work)  and  continuity  within  the  engine  to 


cause 


maintain  equilibrium.  The  engine  rematch  would 

acceleration/deceleration  to  occur  until  equilibrium  was  once  again 
obtained . 
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Ill  Off -Torque  in  Turbofan  Linearization 


Off -Torque 

Off-torque  is  defined  as  any  torque  applied  to  the  high  or  low  spool 
rotor  that  is  different  than  the  torque  required  to  maintain 
equilibrium.  Each  engine  equilibrium  condition  has  associated 
equilibrium  engine  parameters.  The  parameters  of  most  interest  are  Nj  , 

,  T.  ,  W,  ,  and  A.  .  For  the  turbofan  engine,  N,  is  defined  as  the 
low  rotor  speed;  is  the  high  rotor  speed;  T^  is  the  HPT  inlet 
temperature;  W  is  the  fuel  flow  and  A0  is  the  nozzle  exit  area.  The 

I  O 

importance  of  these  particular  parameters  will  be  discussed  shortly. 
Each  of  these  parameters  has  an  equilibrium  value  associated  with  it 
indicated  by  the  subscript  *e".  However,  if  for  example,  was  forced 
to  operate  at  (N^  +  dN^)  and  the  other  four  parameters  remained 
unchanged,  continuity  and  energy  conditions  would  no  longer  be 
satisfied.  The  engine  would  have  a  desire  to  change  both  its  rotor 
speeds  and  HPT  inlet  temperature  to  operate  at  the  new  equilibrium  rotor 
speed  (NJe+dNj).  The  desire  of  the  engine  to  rematch  to  a  new 
equilibrium  point  can  be  used  in  a  computer  simulation  model  cycle  deck 
to  develop  linear  equations  of  motion. 

Off -Torque  Method 

Using  a  cycle  deck  to  force  a  turbofan  to  operate  at  an 

of f -equi 1 ibr lum  point  requires  three  new  iteration  parameters  to  force  a 

new  operating  point  defined  by  (N,  ♦dN,),  N.  ,  and  T,  .  W,  and  A„  are 

le  1  2e  4e  f  8 


not  allowed  to  change.  The  three  iteration  parameters  used  are  low 


speed  spool  off-torque,  Np  high  speed  spool  off-torque,  N2,  and  a 
‘temperature  off-torque',  T^,  made  up  by  adding  dT^  to  T4  to  obtain 

(T  +dT.)  and  multiplied  by  the  turbine  rotor  heat  soak  constant,  r  to 

4  4 

give  T  .  The  heat  soak  constant  is  a  function  of  temperature  and  flow 
and  is  calculated  within  the  cycle  deck.  It  has  units  of  1/sec  and, 
when  multiplied  by  dT^  ,  gives  T^.  The  three  iteration  parameters,  N^, 
N  ,  and  T  are  used  to  define  the  elements  of  the  state  space  equation. 


State  Space  Equation 


The  elements  of  the  A  and  B  matrices  in  the  state  3pace  equation, 

tin  =  [ a] [x]  ♦  [bhuj, 

are  of  the  form. 


i  ,  j  *  1 , 3 


and  B  = 


i= 1 .3;  k= 1 , 2 


where,  for  the  above  turbofan  example. 


X1  =  "l 


U1  =  Wf 


U2  =  A8 


X2  =  N2 


X3  =  T4 


The  off-torques  required  for  operation  at  an  of f -equi 1 ibr lum  point 

are  Mj.  Mj.  and  T^.  Therefore,  the  first  column  of  matrix  A  is  defined 

as  the  off-torques  required  to  rebalance  an  of f-equi 1 ibrium  point  of 

(N.  ♦  dN  )  ,  fL  .  and  T.  for  a  given  W  and  AQ ,  divided  bv  dN, .  The 

le  1  2e  4e  f  8  ‘1 

second  column  of  A  is  similarly  defined  by  the  of f -equi 1 i br lum  point  of 

N.  .  (N„  ♦  d N _ )  ,  and  T  for  a  given  W  and  A  .  The  last  column  of  A 
le  2e  2  4e  °  fe  8 

13  completed  using  (T^  MT^),  and  the  columns  of  the  B  matrix  are 


completed  using  (W  +  dW^)  and  ( Ag  +  dAg)  and  the  off-torques  required  to 
run  to  their  respective  of f-equi 1 ibr lum  points.  Thus,  the  state  space 
equation  linearly  modeling  a  turbofan  engine  about  small  changes  from 
the  equilibrium  is  complete. 

State  Space  Parameters 

It  is  not  initially  evident  why  the  state  space  parameters  Nj  Nj , 
and  are  used  in  describing  the  dynamics  of  a  turbofan  engine.  In 

short,  the  methods  of  choosing  the  proper  engine  parameters  depends  on 
the  speed  of  response  of  the  internal  engine  dynamics,  the  degree  of 
accuracy  required  in  the  state  space  model  and  which  parameters  are  of 
interest.  As  an  example,  pressure  within  the  nozzle,  PQ  becomes 
important  in  describing  engine  dynamics  when  a  long  nozzle  is  involved. 
The  longer  nozzle  allows  the  relatively  high  speed  pressure  dynamics  to 
significantly  affect  transient  conditions.  For  the  F101/B-1B  engine 
used  as  the  example  in  this  paper,  PD  dynamics  are  very  quick  relative 

O 

to  the  length  of  the  nozzle  and  do  not  appreciably  add  to  the  dynamic 
model.  However,  on  some  engines,  like  the  F100/F15 ,F16 ,  P  dynamics  do 

O 

play  a  role  in  modeling  the  dynamics. 

Turbine  inlet  temperature  is  an  important  parameter  in  the  F101/B-1B 
engine  because  of  the  large  HPT  disk  heat  soak  characteristics.  The 
dynamics  of  the  heat  soaked  into  the  disk  significantly  affect 
temperatures  and  pressures  down  stream  of  the  HPT.  Therefore,  the  heat 
soak  time  constant  plays  a  substantial  role  in  engine  transient 
dynamics.  Naturally,  the  parameters  N^  and  N^  are  always  of  interest 
because  the  inertias  involved  with  the  low  and  high  speed  spools 
dramatically  affect  transient  acceleration  times. 


The  size  of  the  deltas  dN^ ,  dN^,  and  dT^  will  vary  in  magnitude  for 
accurate  off-torque  results.  If  the  delta  size  is  too  small  the 
algebraic  limitations  of  the  cycle  would  not  converge  on  an  accurate 
off-torque  value.  If  the  delta  size  is  too  large,  the  off-torque  due  to 
a  positive  delta  would  not  equal  in  magnitude  the  off-torque  due  to  a 
negative  delta.  Delta  sizes  will  vary  for  different  PLA  settings  and 
from  engine  to  engine. 

State  Space  Model  Results 

The  non-linear  transient  engine  cycle  deck  model  for  the  F101/B-1B 
was  used  to  generate  the  state  space  model.  Internal  cycle  deck  logic 
was  developed  to  manipulate  the  Newton-Raphson  convergence  subroutines 
to  use  off-torques  in  rebalancing  to  new  off-equilibrium  points.  State 
space  models  were  developed  at  various  PLA  positions  (equilibrium 
points)  in  small  5  degree  increments  for  accurate  model  continuity 
between  PLA  points.  Figures  3  through  10  show  the  root  locus  and 
gain/phase  plots  for  the  transfer  function  derived  from  the 

various  state  space  models.  The  transfer  function  is  of  most 

interest  in  designing  a  control  scheme  since  the  F101/B-1B  engine  uses 
schedules  to  set  fuel  flow,  ,  for  given  throttle  settings.  Present 
engine  transient  control  logic  uses  as  a  Sain  feedback. 

Notice  the  characteristics  of  the  root  locus  plots.  In  all  cases 
there  are  two  zeros  and  three  poles,  one  pair  as  a  complex  conjugate. 
One  zero  is  relatively  large  and  the  other  situated  near  a  pole.  The 
complex  pole  conjugates  are  located  nearest  the  imaginary  axis.  These 
characteristics  of  relative  root  and  pole  location  are  consistent  at 
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The  task  of  developing  a  control  scheme  could  now  begin,  however,  at 
this  point  in  the  paper,  another  means  of  developing  linear 
thermodynamic  state  space  equations  would  be  of  interest  and  will  later 
be  3hown  to  give  insight  into  the  make  up  of  the  elements  of  the  state 
space  matrices  in  terms  of  temperatures,  pressures,  areas,  air  and  fuel 
f lows . 
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IV  Analytic  Equations 

Analytic  equations  describing  the  thermodynamics  of  a  turbofan 

2 

engine  have  been  derived,  The  equations  were  non-linear,  and  in  all 
cases  failed  to  tie  all  components  of  the  engine  together. 

Recently,  however,  in  a  technical  report  from  General  Electric’s 
Advanced  Engineering  Technology  Department^  an  attempt  was  made  at 
deriving  linear  thermodynamic  relationships  that  encompass  all  stages  of 
a  turbofan  engine.  Although  the  equations  were  not  in  linear 
state  variable  format,  they  did  give  linear  relationships  for  crossing 
various  engine  stages,  the  main  stages  being  the  low  pressure  fan  (LPF)  , 
the  high  pressure  compressor  (HPC) .  the  combustor,  the  high  pressure 
turbine  (HPT)  and  the  low  pressure  turbine  (LPT) . 

To  check  the  accuracies  of  the  linear  analytic  equations,  the  F101 
engine  cycle  deck  model  was  used.  Values  obtained  using  the  linear 
equations  were  compared  to  values  obtained  using  the  non-linear  cycle 
deck  over  small  fuel  step  inputs  corresponding  to  approximately  2 
degrees  PLA. 

Also,  the  equations  derived  below  will  not  be  identical  to  those 
derived  in  the  General  Electric  report.  In  some  cases,  shorter  versions 
of  the  equations  were  found  to  be  less  complicated  yet  described  the 
dynamics  with  greater  accuracy.  It  is  important  to  note  again,  that  the 
partials  derived  by  General  Electric  are  not  in  linear  state  variable 
format.  The  analytic  equations  derived  below  will  keep  linear  state 
variable  format  intact. 


Analytic  Equation  Derivations 

When  describing  the  internal  mass  fiow  rates  and  thermodynamic 


properties  of  a  turbofan  engine,  the  following  equations  describing  mass 


flow  rates,  pressure,  and  temperature  relationships  are  used 
extensively : 

Flow  Function 


where , 


N  = 


wV~T 

AP 


(I) 


W  =  Mass  Flow  Rate  T  =  Temperature  A  =  Area  P  *  Pressure 
N  =  Restriction  Factor  (constant  when  flow  is  choked). 


or  in  terms  of  Mach  number 


N  =  .766  Vy 


y*  l 


2d -y) 


3=  M(2-M) 


(II) 


Pressure  -  Temperature  Relationship 


r- 1 

rn  t 


for  adiabatic  flow  across  the  LPF ,  HPC ,  HPT,  or  LPT  where  y  is  the  ratio 
of  specific  heats  and  r>  is  the  adiabatic  efficiency. 

Static  Pressure  j_Pg)  -  Mach  Number  Relationship 


P 

3 


P  !♦ 


y- 1 


(IV) 


These  equations  link  the  various  engine  components  together  in 
non-linear  relationships.  To  linearize  each  of  these  equations,  the  same 


basic  approach  of  logarithmic  derivatives  is  used. 
For  the  expression 


>^»:i 


>  V- 


A  *  A  *  6k 

9 

B  =  B  ♦  <5B 

c  =  c  +  sc. 

9 

Substitution  gives. 


(A  ♦  6 A)  (B  +  6B) 
e _ e _ 

VC  ♦  6C 

e 


Multiplying  through  and  neglecting  second  order  terms, 

A  B  ♦  A  <5B  +  B  <5A 
e  e  e _ e 

VC~  (1  ♦  6C/C  )  1/2 

e  e 

Substitution  of  the  binomial  expansion 


(1  +  6C/C  ) 


-r—r  £  (1  -  .5(<5C/C  )  )  . 
1/2  e 


and  dividing  by  A  B  rC  gives  the  final  desired  solution  as, 
e  e  e 


6  A  <SB  _o  C 

7T  "  B~  '  -5~ 

e  e  e 


As  stated  before,  the  equations  used  in  this  paper  differ  somewhat 
from  the  equations  described  in  the  General  Electric  report  (different 
intended  uses)  and  will  be  derived  below  in  linear  state  variable  format 
starting  with  the  engine  inlet  and  ending  at  the  engine  nozzle 
algebraicly  connecting  each  engine  stage  to  the  next.  Figure  11 
describes  by  number  the  various  engine  stages  and  locations  of 
pressures,  temperatures,  and  mass  flow  rates  used  in  the  derivation  of 
the  linear  equations  and  should  be  referenced  as  the  equations  progress 
from  3tage  to  stage. 
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Engine  Inlet 

Isentropic  flow  is  assumed  in  the  inlet  giving 


?  -  t 

1  2 


P  =  P 
1  2 


W1  =  W2 


Flow  across  the  fan  is  assumed  adiabatic,  and  the  form  of 
Equation  (III)  states, 


'V1 


r  p 


1  3 


Yint 


1  3 


for  bypass  flow,  and 


'V1 


25 


r  r\ 

c  c 


25 


for  core  f low. 


Taking  the  Log/par tial  and  assuming  a  constant  y  and  efficiency,  T)  in 
the  duct  and  core,  the  equations  become. 


<r<- 1) 

r 


757  K-  ip2)  ■ 


AT13  -  AT 2 


( Bypass ) 


( y  _  -  1 1 

Y 


V’  (AP25-  AP2]  ■ 


c  c 


AT  -  AT 
25  2 


(Core)  . 


If  the  flight  Mach  number  is  held  constant,  AT^  and  AP^  will  be 


r. 

zero.  To 

further 

simplify  the  derivation  of  these 

Lr 

equations  . 

the 

expression  1 Y  -  1 ) /y  •  will  be  defined  as 

5 

GNFD 

f  or 

the 

fan  bypass  flow 

GNF 

f  or 

the 

fan  core  flow 

y 

GNC 

f  or 

the 

compressor  core  flow 

» 

GNT 

for 

the 

turbine  core  flow 
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The  above  equation  for  duct  bypass  flow  and  compressor  flow  are  now 
wri tten  as , 

AP  =  (GNFD)AT  (Duct  Flow) 

1  sJ  i  J 

(2) 

AP  =  (GNF)AT.e  (Compressor  Flow) 

2b  2b 

Fan  Duct  Flow 

The  flow  function  within  the  duct  follows  Equation  (I)  and  gives, 

AN13  =  AWB  ♦  -5A^3  -  A^3  (3a) 

I  f  no  losses  are  assumed  through  the  duct, 

AT13  =  AT16 

AP 1 3  ^  AP  1 6  (3b) 

AN13  =  AN:e 

Compressor  Flow 

Flow  across  the  compressor  is  assumed  adiabatic  satsifying  Equation 
(III)  to  give 

AT3  =  AT25  ♦  AP3/GNC  -  AP^/GNC  (4) 

Engine  Combuster 

The  equation  governing  the  flow  and  heat  gain  in  the  combustor  is 

W,_  (h,-  -  h,)  =  W. (LHV)  (5) 

2b  oy  o  t 

where , 


25 

h 

Wf 

LHV 


1  39 


air  flow  in  core 

enthapy  at  specified  engine  iocation  (equal  to  Cp»T) 
fuel  flow 

lower  heating  value  of  fuel  (assumed  constant) 

turbine  inlet  temperature  prior  to  turbine  rotor  heat  soak 
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The  turbine  rotor  on  the  F101  engine  13  very  massive  and  the 


temperature  heat  soak  dynamics  play  a  sizable  part  in  determining 
engine  dynamcs . 

Taking  the  log/der lvative  gives 

AW.  =  AW„C  +  AT,. (Z . / (Z  - 1 ) )  -  AT,/(Z,-1)  (5a) 

f  25  39  4  4  34 

where  Z 4  =  T^/T^ 

Across  the  combustor,  the  pressure  drop  is  assumed  small  enough  and 
the  pressure  ratio  P^/P^  nearly  constant  such  that, 

AP3  =  (5b) 

The  accuracy  of  this  assumption  will  be  discussed  in  the  error  analysis 
section . 

Turbine  Inlet 

Flow  into  the  turbine  inlet  is  assumed  choked.  This  may  not. 
however,  be  the  case  at  near  idle  conditons.  but  using  the  cycle  deck 
model,  the  restriction  factors  could  be  computed  at  the  low  idle  speeds. 
For  the  majority  of  engine  operations, 


- - - - =  Constant 

#41P41 


AP,  ,  AW.C  *  .  5  AT , .  (6) 

4  25  4 

where  AT4  has  been  previously  defined  in  the  empirical  equation 
discussion  as  TIT.  TIT  and  T^  can  be  used  interchangeably. 

High  Pressure  Turbine 

The  assumption  made  for  flow  across  the  high  pressure  turbine  is 

38 


v.V. 


>  .V. 


%  -N 


This  is  a  reasonably  accurate  assumption  and  will  be  discussed  later 


during  the  error  analysis  section. 

Low  Pressure  Turbine 

As  before,  adiabatic  flow  across  the  low  pressure  turbine 
assumed.  The  appropriate  equation  is. 


which  gives, 

APg  =  Ap42  ♦  (A^  -  A]j  )/GNT. 

The  equation  for  the  restriction  factor, 


W  =  -/r 
5  i 

a5p5 


gives  , 


(7a) 


an5  =  aw25+  ,5AT5  -  ap5 


(7b) 


Unbalanced  Power 


Fan 

The  horsepower  oi  the  fan,  (Fp) ,  is  given  by  definition  as. 
FP  =  (T.,  -  T . ) W  , 

i  sj  2  2 

or  after  the  log/der lvative  . 


is 


AFP 


AW.  *  AT  Z. 


(8) 


where,  Zg  =  Ti3/ <Tl3_T2 > ' 

Low  Pressure  Turbine  (LPT) 

The  horsepower  of  the  low  pressure  turbine,  (LP) ,  is  given  by 
definition  as, 

LP  =  (T.  -  T  .„) W  ,  which  produces, 

LP  *  AW25  ♦  AT42Z9  ♦  ATg(l-Z9)  (9) 

where,  Zg  =  T42/(T42~T5) 

The  difference  between  the  fan  horsepower  and  the  LPT  horsepower 
would  be  zero  in  equilibrium.  However,  when  not  in  equilibrium,  the 
differences  are  equal  to  the  off-torque,  I^Nj,  where  I ^ i s  the  low 
pressure  spool  inertia.  In  log/derivative  format, 

AN  =  ALP  -  AFP 

'  AW25  +  AT4Z9  +  AT5(1“Z9)  -  AW2  +  AT13‘  (10) 


Compressor 


In  a  similar  manner  to  the  LPF  and  the  LPT,  the  compressor 
horsepower,  CP,  log/differential  is, 


acp  -  aw25  ♦  at3z7  ♦  at25  (1*Z7) 

where.  Z?  =  T3/(T3-T25). 


(11) 


High  Pressure  Turbine  (HPT) 


The  equations  governing  the  HPT  are 


V.V.  V.  AA'f.  ^  r»''LTV.'Tk.’VV  V  V 


HP  =  (T42  -  T4)W25  and' 


AHP  =  AW  ♦  AT 

25  4 


remembering  that  AT^  *  AT^ . 


The  difference  in  the  torque  for  the  HPT  and  compressor  gives  the 
off-torque  for  the  high  speed  rotor  as, 


AHP  -  ACP  =  AT  -  AT,Z_  -  AT  C(l-Z_). 

4  37  25  7 


Using  Equations  (2),  (4),  (5b),  and  (6),  the  AT^  term  can  be 
rewritten  as, 

AT,  =  AT,_ ( 1  -  GNFD/GNC)  ♦  W„,/GNC  ♦  TJ/(2(GNC))  (14) 

3  25  25  4 

Substitution  of  Equation  (14)  into  Equation  (13)  gives, 

AN2  =  AHP  -  ACP 

=  -  AW  Z„/GNC  +  AT  ( 1  -  Z_/ ( 2 (GNC) )  +  AT  (Z  (GNFD/GNC)  -1). 

25  7  4  J  25  / 


Temperature  Heat  Soak 


In  equilibrium.  T^g  is  equal  to  T^  since  the  HPT  rotor  is  in  thermal 

equilibrium.  When  fuel  is  added.  T,  immediately  rises,  However,  T 

oy  *1 

takes  some  finite  amount  of  time  to  regain  thermal  equilibrium.  The 
differential  equation  which  describes  the  turbine  rotor  thermal  dynamics 
due  to  fuel  flow  input  is 


T4  =  K4Z3(T39> 


where  T,n  has  already  been  shown  a  function  of  W„_ ,  T  ,  T„_  and  W, 

3SJ  25  4  25  f 

Taking  the  log/deri vative ,  the  equation  becomes, 


.\V. 

WvV  W  "> ' 

*  O  V  •  •  O'V I  V 


(16) 


AT4  =  VSV 


where  is  the  heat  soak  time  constant  and  Z^T^/T^.  Using  the 


expression  for  ATjg  from  Equation  (5a)  and  substitution  of  Equation  (14) 


into  Equation  (16)  gives, 


AT  .  /K  =  (Z./Z.)  [AW„C  (  ( 1  -ZJ  +  l/GNC)  ♦  AT  .  /  (2  (GNC) )  + 

1  4  4  4  25  4  4 


AT, ( 1 -GNFD/GNC)  ♦  AW  (Z.  -  1 ) ] 
a  D  2  4 


Equations  (10),  (15),  and  (17)  make  up  a  partial  state  space  matrix 


of  the  form, 


[X]  =  [A] [W]  +  [B ] [U]  ♦  [ E ] [ V ] 


or  in  terms  of  engine  parameters, 


AT  /V 

4  4 


[A]  AW  ♦  [B]  f  +  [ E ]  AT 


where , 


[A]  = 


-z?/gnc 


( 1 -Z?) / (2 (GNC) 


Z,/Z. ( 1 /GNC ♦ (1-Z.) 

3  4  4 


1/(2 (GNC )  ) 


0  (1-Z, 


[  B  1  =  I  0 


0  .  [E]  =  0 


Z7 (GNFD/GNC) -  1 


(Z.-l)  0 

4 


0  0 


( 1 -GNFD/GNC) 


To  transform  the  above  equation  into  the  conventional  form, 


txi  =  mm  +  cb  j  cu] . 

the  continuity  equations  for  a  turbofan  engine  are  employed. 

The  total  nozzle  area  for  the  F 1 0 1  is  made  up  of  the  bypass  duct 


area,  A„n,  and  the  core  area  A  .  Remembering  the  restriction  factor  in 
8B  oC 


the  duct,  N  is 

N 


13  *13Pi3 


(19) 


continuity  of  flow  at  station  8  (the  nozzle)  gives. 


N  1  3  A  1  3  P  1  3 


VT 


1  3 


W8A8CP8 


Taking  the  log/der i vati ve  and  recalling  AA13  is  zero, 

AA8B  =  AN13  +  '  5AT0  -  -5AT13  '  ^8  +  ^13-  <20) 

Flow  at  station  5,  the  LPT  outlet,  must  also  satisfy  continuity. 

The  relationship  gives, 


AA8C  =  AN5  +  -5AT8  -  -5AT5  '  ^8  +  AP5-  (21) 

As  stated  earlier,  the  nozzle  area  is  made  up  of  both  core  and 

bypass  area  such  that  A  =  A  +  k  ,  or  after  the  log/derivati ve , 

O  OD  Ot 

A8AA8  =  *88  A^3B  +  '  (22a) 

Each  area  satisfies  the  following. 
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(22c) 


A 


8B 


W 


il 

P 


A- 


( 2  2d ) 


Dividing  Equation  (22b)  by  (22d)  gives, 

AAg  ( 1  +  B )  =  Ai^gB  ♦  A .  (23) 

where  B  is  the  air  flow  bypass  ratio  W,_/Wn_.  Substitution  of  Equations 

1 0  ^0 

(20)  and  (21)  into  Equation  (23)  gives, 

AA  (1  +  B)  =  AIL  B  *  .5A^B  -  .5A^  B  -  APg  B  ♦  AP^B  + 

AN5  +  .  5ATq  -  ,5A^  -  A|  ♦  A^  .  (24) 

When  two  seperate  air  flows  mix,  in  this  case  the  duct  bypass  and 
core  air  flows,  an  approximation  of  the  final  temperature  state  is  a 
flow  weighting  of  the  two  air  flows.  At  station  8,  the  mixing  point, 
the  equation  for  temperature  is, 


(W 


25 


W13)T8  S  W13T13  +  W25V 


Rewriting  in  terms  of  constant  bypass  ratio,  B,  and  taking  the 
log/derivative  . 


AV‘  *  91  '  AT13  *  A15  ' 

where  Q  =  BlTj^/T^. 


(25a) 


A  simplified  means  of  mixing  pressures  P^3  and  P^  based  on 
experimental  results  of  tested  turbofan  engine  data  gives  the 


log/den vative  relationship, 


(25b) 


APq  ( 1  +  B)  =  AP^B  ♦  APg. 

Accuracies  of  the  assumptions  made  in  Equations  (25a)  and  (25b)  will 
be  considered  in  the  discusion  on  error  analysis.  Substituting  Equations 
(25a)  and  (25b)  into  Equation  (24)  results  in, 

AA_  ( 1  +  B)  =  All  B  +  AN  ♦  .5AT.Z.  -  .5AT..Z..  (26) 

8  Id  5  05  lo  5 

where  Zc  =  (B  —  Q) / ( 1  +  Q) . 

Substitution  of  Equations  (7a)  and  (7b)  into  Equation  (6)  defines 
the  restriction  factor  at  station  5  as 

AN5  =  AT4(1/GNT  -1/2)  -  A^  (1/GNT  -1/2).  (27 

Further  substitution  of  Equations  (2),  (3b).  and  (27)  into  Equation  (26) 
gives  the  first  continuity  equation  in  final  form  being, 

(1  ♦  B)  -  A^5  ♦  AT4  (1/GNT  -  1/2)  + 

AT12(B/2  -  B(GNF)  -  Zg/2)  ♦  A1^  (^ /2  -  1/GNT  ♦  1/2) 

(28) 


AA8  ( 1  ♦  B)  =  A^ 


Another  continuity  conditon  must  also  be  satisfied.  When  the  duct 
and  core  mass  air  flows  mix,  the  static  prssures  at  the  mixing  point 
must  be  equal,  or  for  the  F101  engine. 


P  =  P 
s  16  s5 

Recalling  that  there  are  no  looses.  P  ,,  is  also  ecual  to  P  just  as 

sib  s  1  o 

P  is  equal  to  P,,.  These  terms  will  be  used  interchangable . 

The  relationship  between  restriction  factor  and  Mach  number,  M,  is 
given  by, 


N.  . 

1 6 


'  6  6  V v  M 


16 


(  29 ) 


This  equation  becomes  very  cumbersome  when  taking  the  1 og / der l vat l ve . 


Based  on  past  General  Electric  performance  data  from  turbofan  engine 
testing,  a  simplified  equaton  for  the  restriction  factor  over  the 
operational  range  of  an  engine  is, 

»,6  =  V2  •  Ml6>  <38* 

The  accuracy  of  this  equation  relative  to  Equation  (29)  can  be  seen  in 
Figure  12  which  plots  Equation  (30)  versus  Equation  (29)  as  a  function 
of  Mach  number.  The  relationship  between  static  and  total  pressures  as 
a  function  of  Mach  number  is  by  definition, 


16 


16 


s  16 


r , 


i  + 


16 


M 


16 


^16"  1 


(31a) 


Taking  the  log  of  both  sides,  and  using  the  log  senes  approximation, 
Equation  (31a)  becomes, 

*°®Ps 16  =  Io«PI6  •  -5p16M?6  ,3Ibl 

Considering  station  16  and  taking  the  log/der i vative ,  Equations 
( 30) ,  and  (31b)  yield , 


iNi6  ’  “-V211  •  “l6,/l2  -  V  ,32»> 

ap=i6  1  APie  -  l32b) 

Combining  Equations  (32a)  and  (32b)  gives  the  desired  result  in  terms  of 
static  pressure  as 

APS16  =  Apn  -  Af^  ,32c) 

2 

where  P  is  substituted  for  and  2.  =  Y M../2(2  -  M.c)/(1  -  M,c). 

10  10  110  10  1 0 

Similarly  for  station  5. 


APs5  =  AP5  -  \ 


( 3  2d ) 


46 


where  Z,  =  .5M,(2  -  M_ ) / ( 1  -  M.)  . 

2  5  5  5 

Since  static  pressures,  AP^  and  APg^3  are  equal, 

^13  -  AN16Z1  =  ^5  -  AN5Z2-  (33) 

Recalling  the  relationships  in  Equations  (1),  (2),  (6)  and  (7a), 

Equation  (33) ,  can  be  written  as 


AT,, ( ( 1+Z, )GNFD  -  Z,/2)  -AT, (Z./GNT- . 5)  ♦  1/GNT) 

1  o  I  1  5  2 

=  AW.Z, ( 1 +B) /B  ♦  AW  ( 1  -  1 / B )  -  AT. ( 1/GNT- , 5)  (34) 

2  i  25  4 

This  is  the  final  form  to  be  used  in  the  second  continuity  relationship. 

Finally,  the  relationship  between  T  and  Tn,  can  be  determined 

lo  25 

experimentally  using  the  engine  cycle  deck  at  the  desired  flight 
condition.  For  a  given  flight  condition. 


where  power  factor,  Kf  ,  is  the  factor  relating  T13  and  T25  derived  from 
the  cycle  deck.  Taking  the  log/der iviati ve , 

AT13  =  AT25(K{).  (35) 

Equations  (28),  (34),  and  (35),  when  combined,  give  an  important 
transformation  matrix  of  the  form, 

C  G  ]  C  V  ]  =  CH]  CW]  ♦  [JHU],  (36) 

where . 


Rewriting  Equation  (36)  as, 

[V]  =  [Gi ~ 1 [H H W J  ♦  CG ] " 1 C J ] [ U ]  (38a) 

Equation  (38a)  will  be  substituted  into  Equation  (18)  in  place  of 


the  [V]  matrix. 


The  other  transformation  needed  is  between  the  [X]  and  [W]  matrix. 


This  is  obtained  from  the  fan  and  compressor  maps.  Fan  and  compressor 


maps  are  experimentally  derived  functions  of  flow  versus  pressure  ratio 


along  constant  rotor  speed  lines.  The  cycle  deck  uses  these  map3  to 


define  an  operating  point  (an  air  flow  and  pressure  ratio).  From  these 


maps,  around  small  equilibrium  operating  points,  air  flow  can  be 


considered  a  linear  function  of  pressure  ratio  and  speed.  For  the  fan 


and  compressor  this  means, 


AW2  =  AN^  (DW2QN)  ♦  AP  (DWQP) 


AW25  =  A^5  (DWCQN)  ♦  AP3QP2g(DWCQP) 


where  DW2QN.  DWQP,  DWCQN,  and  DWCQP  and  linear  derivatives  derived  from 


the  cycle  deck.  From  Equation  (2),  (5b),  (6),  and  the  definition  that 


AP  QP  =  AP  -  AP 
3vr25  3  25' 


the  equations  can  be  rewritten  as 


AW2  =  AN  (  D  W2Q  N )  +  ATJ;J  (DWQP)  (GNFD) 


AW25  *  pWggp  AN2  (DWCQN)  ♦  AT4  (DWCQP) 


-  AT  _ (DWCQP) (GNF) 
/o 


In  matrix  form,  the  above  equations  are, 


AW25  =  [ DWX]  AN2 


♦  [ DWV ]  AT, 


(38b) 


AT./K. 

4  4 


which  are  of  the  form. 


[W]  =  [DWX] [X]  *  [ DWV ] [ V ] 


where , 


\a> 

vv 

v*V 

s-\'- 
v  v 


5?? 


v.v.t 
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CTO 


CTO 


vv;.-V-> 


r  DW2QN 


0 
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DWX  = 


DWCQN 


( 1 -DWCQP) 


K  (DWCQP) 

4 

2 ( 1 -DWCQP) 
K. 


and 


DWV  = 


DWQP(GNFD)  0 
0  0 

0  0 


1 


( 1 -DWCQP) 


[ -DWCQP (GNF) ] 


The  transformation  matrices  of  Equations  (38a)  and  (38b)  can  now  be 
substituted  into  Equation  (18)  to  obtain  the  conventional  form  of, 

[X]  »  [ A] [X]  ♦  [ B ] [ U ]  . 

The  final  form  of  Equation  (18)  after  the  transformation  substitutions 
is  , 

[X]  *  (A  ♦  EG" 1 H) C I  -  [DWV] (G_1H] ] [DWX] [X]  ♦ 

- !  - !  - !  - !  (39) 

[(A  ♦  EG  H)  [  I  -  [DWV]  (G  H)  ]  [DWVHG  J)  ♦ 

(B  +  EG'^)  ][U]  , 

where  [I]  13  the  identity  matrix. 

Equation  (39)  is  in  the  proper  state  space  format. 


Error  Analysis 

The  error  analysis  will  include  consideration  of  the  assumptions  in 
the  development  of  the  analytic  state  space  model  and  errors  induced 
into  the  model  with  the  use  of  the  transformation  matrices. 

One  important  assumption  used  in  the  development  of  the  analytic 
equations  was  that  AP  is  equal  to  AP  This  assumption  is  based  on 


the  ratio  P  /P  being  nearly  constant  over  the  range  of  turbofan  engine 

operation  of  concern.  When  the  log/derivative  is  taken  the  assumed 

result  is  achieved.  The  other  important  assumption  made  was  that  AT^ 

the  HPT  inlet  temperature,  is  equal  to  AT^.  inlet  temperature. 

This  assumption  is  based  on  the  nearly  constant  HPT  efficiency  over  a 

wide  range  of  engine  operating  conditions.  The  assumption  was  tested 

using  engine  cycle  deck  data.  Over  the  range  of  35  degrees  PLA  to  70 

degrees  PLA,  the  ratio  of  AT.  to  AT.„  was  .983  ±  .015  verifing  the 

4  42 

accuracy  of  the  assumption. 

The  assumptions  made  in  deriving  Equations  (25a)  and  (25b)  were 

based  on  past  engine  experience  and  common  practice  when  combining  two 

separate  air  flows  like  W^  and  W^.  The  accuracy  of  mass  weighting 

temperatures  when  combining  air  flows  is  seen  in  Figure  13a  which  plots 

T  calculated  using  the  cycle  deck  versus  T  approximated  by  mass 
o  o 

weighting.  The  slope  of  the  line  is  nearly  one  giving  strength  to  the 
approximation.  However,  when  the  log/derivative  is  taken,  the  accuracy 
of  the  resulting  Equation  (25a)  is  reduced  to  approximately  80%  over  the 
range  of  engine  operation.  Equation  (25b),  however,  proves  much  more 
accurate  in  its  approximation  of  pressure  addition  of  two  flows.  Figure 
13b  shows  a  plot  of  AP  as  calculated  by  the  cycle  deck  versus  AP  as 

O  O 

calculated  by  Equation  25b.  Again,  the  slope  is  nearly  one.  confirming 
the  accuracy  of  the  approximation. 

The  largest  error  in  computing  the  elements  of  the  analytic  state 
space  model  come  from  the  numerous  substitutions  and  transformation 
matrices  used  to  reduce  the  analytic  model  to  a  3  by  3  state  space.  The 
errors  are  large  enough  to  require  the  user  to  rely  on  the  empirically 
derived  state  space  models  for  accurate  modeling  of  engine  transient 


DELP8  DCK 


characteristics  and  the  analytic  models  for  insight  into  individual 
engine  parameter  effects  on  transient  performance 

Comparison  of  the  analytic  versus  the  empirical  state  space  models 
show  the  substantial  differences.  Keep  in  mind,  however,  the  form  of 
the  empirical  equations  and  the  form  of  the  analytic  equations.  The 
states  of  the  analytic  equations  were  defined  as 

M.if  . 

e 

For  comparison,  it  is  necessary  to  multiply  row  1  of  the  [B]  matrix 

of  Equation  (39)  by  N.  ,  row  2  by  N„  and  row  3  by  T.  .  Then, 

le  /e  4e 

multiplying  columns  1  and  2  of  the  [B]  matrix  by  W,  and  A„ 

fe  8e 

respectively,  will  put  the  analytic  form  of  the  equations  into  the 

empirical  form. 
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V  Engine  Control  Design 

Using  the  analytic  equations  described  above,  state  space  equations 
were  derived  at  various  PLA  settings  starting  with  35  degrees  to  70 
degrees  in  5  degree  increments.  Table  I  compares  the  analytic  state 
space  models  with  the  previously  derived  empirical  state  space  models 
for  each  PLA  setting.  The  errors  induced  by  numerous  substitutions  and 
transformation  matrices  do  not  give  good  correlation  between  to  two. 

There  is  one  important  insight  the  analytic  state  space  models  do 

give.  The  [B]  matrix  of  the  analytic  equations  shows  that  T^  is  the 

only  parameter  affected  directly  by  fuel  flow  inputs.  The  empirical 
state  space  models,  however,  show  small  effects  of  fuel  flow  input3  on 
Nj  and  N^.  Upon  closer  examination  of  the  root  locus  plots  generated 
from  the  empirical  equations,  (see  Figures  3  through  10),  a  large 
negative  zero  is  generated.  If,  in  the  empirical  equations,  the 
elements  B(l,l)  and  B ( 1 , 2 )  of  the  [B]  matrix  are  set  equal  to  zero 
corresponding  to  analytic  values,  the  large  zero  characteristic  in  each 
of  the  the  root  locus  plots  disappears.  The  loss  of  the  large  zero  does 
not  affect  the  character istcs  of  the  empirical  state  space  models. 
However,  due  to  the  differences  in  the  empirical  and  analytic  values,  it 

is  necessary  to  use  the  empirical  form  in  developing  feedback  control 

logic.  The  analytic  equations  can  still  be  used  to  understand  trends 
and  the  effects  of  specific  parameters  on  engine  performance. 

The  reason  for  the  large  negative  zero  is  understood  when  examining 
the  accuracy  of  convergence  routines  used  today  in  non-linear  engine 
cycle  decks  to  balance  energy  and  continuity  as  discussed  previously. 
The  convergence  subroutine  works  by  equating  changes  in  N ^ ,  T  ,  W 
and  Ag  with  off-torques.  and  N^,  and  T^.  The  convergence  subroutine 
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does  not  know  that  and  are  not  directly  affected  by  W^  and,  by  its 
nature  must  generate  small  off-torque  values  and  N^,  for  a  given 

change  in  W  .  Although  the  error  is  small  it  manifests  itself  in  a  large 
negative,  non-effective  root. 

Using  the  empirically  derived  state  space  matrices  between  35 
degrees  PLA  and  70  degrees  PLA,  feedback  control  systems  were  derived  to 
improve  the  F101  engine  response  to  small  changes  in  PLA. 

The  present  control  logic  internal  to  the  engine  cycle  deck 
simulation  involves  the  feedback  of  the  fuel  metering  valve  rate,  and 
N  .  For  PLA  changes  of  2  degrees,  response  is  on  the  order  of  3  seconds 
to  reach  steady  state  values. 

The  present  feedback  control  logic  installed  in  the  engine  cycle 
deck  had  to  be  bypassed  and  a  new  feedback  logic  using  a  lead/lag  filter 
equivalent  to  a  N2  and  N2  feedback  was  installed.  The  lead/lag  filter 
was  used  to  control  to  a  desired  gain  margin  for  tracking  error  and  to 
increae  phase  margin  for  increased  engine  response  time.  With  a  new 
control  feedback  system  based  on  the  empirical  state  space  models, 
response  was  improved  up  to  75%  in  rise  time.  Figure  14  illustrates  the 

new  feedback  system.  N was  used  specifically  for  the  F101  because  the 

F101  is  controlled  by  a  N2  control  schedule  which  is  specified  by  a  PLA 
setting.  In  an  engine  like  the  F 1 1 0 ,  which  is  controlled  by  an 
schedule  based  on  a  PLA  setting,  N.  and  could  have  been  fed  back. 

Table  II  tabulates  the  gains  on  the  rate,  ,  and  Nj  feedback  for 
various  PLA  settings.  Also  listed  in  Table  II  is  the  percent  increase 

in  response  times  to  90%  steady  state  values.  Figures  15a  through  1 5 h 

show  the  increases  in  the  response  of  over  the  present  transient 


cvcle  deck  model . 
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IV  Conclusions 

The  analytic  equations  offer  a  means  of  visualizing  the  effects  of 
individual  parameters  on  engine  performance  but  are  too  numerically 
inaccurate  to  be  used  in  modeling  transient  engine  performance. 
Emperical  equations,  however,  offer  accurate  state  space  models  for  the 
developement  of  engine  feedback  control  logic. 

Using  the  empirical  state  space  models  to  develop  new  feedback 
control  logic,  the  F 1 0 1  transient  performance  was  improved  over  present 
engine  performance  for  small  changes  in  PLA. 
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